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FOREWORD 

This  preliminary  report  was  written  to  provide  a reference  to  anyone  interested  in  utilizing  the 
optical  design  and  analysis  computer  program  written  by  W.  C.  Fitzgerald.  The  writing  of  this  report 
was  done  in  June  and  July  of  1975  and  it  was  reviewed  for  technical  accuracy  by  W.  C.  Fitzgerald. 
This  program  has  been  used  in  the  design  and  analysis  of  several  systems  including  AIM-9L,  AGILE, 
and  APD  seeker  systems. 

This  report  is  released  at  the  working  level  for  information  only. 


P.  G.  ARNOLD,  Head 
Electro-Optics  Division 
Weapons  Department 
25  February  1976 
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INTRODUCTION 


Manually  performing  exact  numerical  ray  traces  for  complicated 
optical  systems  Is  always  a tedious  task.  Accordingly,  simplifications 
' and  approximations  can  be  applied  to  arrive  at  paraxial  ray-tracing 
techniques.  These  techniques  are  quicker,  but  are  frequently  Inadequate 
- to  describe  real  optical  systems. 

However,  modern  digital  computers  can  be  programmed  to  conduct 
exact  numerical  ray  traces  and  even  to  optimize  system  parameters  when 
a fairly  close  approximation  to  the  ideal  system  is  read  in  as  data. 

Such  a program  is  342*WCRF. OPTIC,  designed  by  W.  C.  Fitzgerald  of  the 
Naval  Weapons  Center  (NWC) . This  report  describes  the  program  and  gives 
Instructions  for  using  it. 


The  user  knowledgeable  in  optics  can  apply  this  program  to  a large 
variety  of  optical  and  Infrared  design  problems,  including  the  determina-  ' 

tion  of  static  gain  curves  for  reticle-modulated  seeker  systems. 

The  purpose  of  this  report  is  to  help  the  optical  designer  with 
little  or  no  knowledge  of  programming  to  make  efficient  use  of  this 
powerful  designing  tool.  To  this  end,  a flow  chart  (Figure  1) , an 
outline  of  capabilities  and  options  open  to  the  user,  and  Instructions 
for  implementation  of  the  options  are  included.  In  addition,  definitions 
of  the  aberrations  used  in  the  program  and  a discussion  of  various  methods 
of  calculation  appear  within  the  text.  A sample  problem  is  run  and  the 
program  output  is  interpreted,  and  finally  a troubleshooting  section 
explaining  some  common  error  statements  is  Included.  . 4 
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CAPABILITIES  AND  OPTIONS  OPEN  TO  THE  USER 


PROGRAMMING  OPTIONS  OF  342*WCRF.OPTIC 

1.  The  user  can  choose  batch  or  demand  mode  of  operation. 

2.  Input  data  is  used  to  select  various  computation  options, 

as  follows : 

a.  Spherical,  aspheric,  aconic,  parabolic,  hyperbolic, 
and  elliptical  surfaces 

b.  Refracting  or  reflecting  surfaces 

c.  Tilts,  decenterings 

d.  Ray-trace  only  mode,  design  mode,  static  gain  mode 

e.  Graphics  options:  Ray  trace  or  spot  diagram  (other 
versions  of  WCRF .OPTIC  may  be  required  for  CRT  display 
of  graphics) 

f.  Whether  ray  traces  for  1,  2,  or  3 colors  are  desired 

g.  Numbers  of  and  parameters  of  rays  to  be  traced  through 
system. 

3.  Input  data  can  also  control  program  flow  of  operations , as 

follows : 

a.  Several  operations  on  Input  data  may  be  performed  within 
one  run. 

b.  Operations  on  the  results  of  previous  calculations  may 
be  performed  within  the  same  run. 

c.  Several  different  problems  may  be  computed  within  one 
run. 


RAY-TRACE  MODE 

Used  to  demonstrate  performance  of  a given  system. 

x.  Any  numbers  of  rays  (up  to  400  for  1 color,  900  for  3 colors) 
with  any  specifications  may  be  traced  through  as  many  as  30 
Interfaces. 
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2.  Fans  of  rays  or  rings  of  rays  can  be  generated  (user  specifies 
number  of  rings  or  number  of  rays  in  fan) . 

3.  Graphics  options  Include: 

a.  No  plots 

b.  Ray  tracing  (user  specifies  graph  axis  extremes) 

(1)  Plane  for  traces,  X-Y,  X-Z,  Y-Z 

(2)  Which  surface  of  the  system  the  diagram  starts 
with 

(3)  New  or  old  plotting  frame 

c.  Spot  diagrams  (user  specifies  graph  axis  extremes) 

(1)  Plane  for  spot  diagram  X-Y,  X-Z,  Y-Z 

(2)  Interface  nunber  of  spot  diagram 

(3)  Option  to  find  image  plane  and  image  size  and 
make  spot  diagram  at  this  location 

(4)  New  or  older  plotting  frame. 

4.  Data  can  be  written  for  a specified  surface  on  to  the  last 
surface . 

5.  The  nunber  of  rays  within  a specified  Incremental  radius 
about  the  focused  spot  may  be  requested. 


DESIGN  MODE 

(Used  to  optimize  a given  system  by  varying  system  parameters)  . 

1.  User  may  specify  N variables  which  the  program  may  adjust 
(radii,  position,  etc.) 

2.  User  then  selects  up  to  N of  13  aberrations  to  be  minimized. 

3.  Relative  weighting  may  be  applied  to  aberrations  desired  to 
be  corrected. 

4.  Maximum  nunber  of  iterations  may  be  specified. 

3.  Limits  on  increments  arc  determined  by  the  user. 
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6.  Output  results  Include: 

a ■ Aberration  terms  and  variable  parameter  values  corres- 
ponding to  Iteration 

b.  Ray  data  for  17  rays  traced  through  the  system  defined 
by  last  lnteratlon 

c.  Graphical  ray  trace  of  system  (spot  diagrams  not 
possible  in  the  design  mode) . Options  are  the  same 
as  ray-trace  mode  graphics  options  for  ray-tracing 
plots. 


STATIC-GAIN  MODE 

In  this  mode  information  generated  by  the  calculation  of  a spot 
diagram  is  used  to  predict  the  modulation  efficiency  versus  field  angle 
for  a rotating  mirror/reticle  system. 

1.  Since  spot  diagram  information  is  required,  the  ray— trace 
mode  must  be  used  in  conjunction  with  the  static-gain  mode. 

i*  2.  Any  reticle  configuration  can  be  read  into  the  system  as  data. 

3.  Any  frequency  response  curve  can  be  read  into  the  system  as 

. data. 
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INSTRUCTIONS  FOR  IMPLEMENTING  SPECIFIC  MODES 


OF  OPERATION  AND  OPTIONS 


GENERAL  INFORMATION 

Calling  the  program  is  accomplished  by  the  following  Instructions: 


@RUN  (Proper  user  identification) 

0ASG.A  342*WCRF. 

@XQT  342*WCF. OPTIC 
. . . .DATA. ... 

@FIN 

The  data  whether  on  file  or  cards  will  be  of  the  following  general 

form : 

1.  User  name  and  Job  name  card  (for  graphics  output  identification) 

2.  System  and  graphics  specification  card 

3.  Surface  data  and  ray  data  cards 

4.  Type  of  computation  to  be  conducted  cards. 

To  continue  with  another  problem  within  the  same  run  it  is  neces- 
sary to  Include  the  following  after  the  type  4 cards  above: 

1.  System  and  graphics  specification  cards 

2.  Surface  and/or  ray  data  cards  (only  for  parameters  which  the 
user  wishes  to  change) 

3.  Type  of  computation  to  be  conducted  cards. 


The  conventions  adopted  are  as  follows: 

1.  Mirrors  shall  have  Indexes  of  refraction  of  -1.0  at  all  colors; 
Mangln  mirrors  shall  have  indexes  of  refraction  of  -n. 

2.  The  X-axis  corresponds  to  the  optic  axis,  Y is  up,  and  Z Is 
out  of  the  page. 

3.  A parabolic  or  aconic  surface  is  denoted  by  a normal  surface 
card  followed  by  a surface  card  containing  the  negative  of 
the  surface  msnber  and  the  aspheric  constants. 
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4.  Instead  of  using  on  axis  lens  thickness,  this  information  is 
contained  in  the  input  in  the  form  of  the  location  of  the 
surface  along  the  X-axis. 

5.  Distances  may  be  expressed  in  any  unit. 

6.  Radii  of  curvature  are  positive  when  the  center  of  curvature 
occurs  on  the  right  side  of  the  surface. 

7.  Elliptic  and  hyperbolic  surfaces  are  denoted  by  a surface 
card  with  surface  number  greater  than  100  plus  the  necessary 
defining  data. 


DATA  CARDS 

Lists  of  items  to  be  included  on  data  cards  for  formats  and  option 
specification  are  given  below.  Refer  to  the  coding  sheet  in  Figure  2 for 
column  information  and  examples. 

Plotting  Identification  Card 

Plotting  identification  card  must  include  NWC  code  number  and 
user  name  and  can  Include  job  name.  Up  to  36  alphameric  characters  can 
be  used.  Refer  to  Example  1 on  FORTRAN  coding  sheet  for  format 
(Figure  2). 


System  and  Graphics  Specification  Caid 

R®fer  to  Example  2 on  FORTRAN  coding  sheet  for  format  (Figure  2) . 

1.  IT  “ Nunber  of  Interfaces  (surfaces) 

2.  JT  ■ Number  of  rays 

3.  IW  ■ First  interface  for  which  data  is  to  be  written 

4.  NR  * Specifies  type  of  ray  diagram  desired 

a.  NR  - 1:  No  ray  diagram  (NR  must  • 1 if  NS  j*  1) 

b.  NR  - 2:  X versus.  Y plot 

c.  NR  - 3:  X versus  Z plot 

d.  NR  » 4:  Y versus  Z plot 

5.  1R  ■ First  Interface  to  appear  on  ray  diagram  • 
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6.  NS  - Specifies  type  of  spot  diagram  desired 

a.  NS  ■ 1:  No  spot  diagram 

b.  NS  - 2:  X versus  Y diagram 

c.  NS  «*  3:  X versus  Z diagram 

d.  NS  ■ 4:  Y versus  Z diagram 

e.  IS  - Interface  number  of  spot  diagram 

f.  NF  ■ 1 -*•  New  plotting  frame;  = 2 -*•  same  plotting  frame 

g.  XL  = Lower  limit  of  x values  to  appear  on  plot 

h.  XU  ■>  Upper  limit  of  x values  to  appear  on  plot 

i.  YL  = Lower  limit  of  y values  to  appear  on  plot 

J . YU  ■ Upper  limit  of  y values  to  appear  on  plot 

NOTE:  If  Y-Z  plane  is  requested  for  either  ray  diagrams  or  spot 
diagrams,  XL,  XU  refer  to  lower  and  upper  y values  and  YL,  YU  refer  to 
lower  and  upper  z values. 


Surface  Data  Cards 

Refer  to  Example  3 on  FORTRAN  coding  sheet  (Figure  2) . 

N1  - 1 Read  interface  data 

a.  0 < N2  < 100  spherical  surface 
N2  ■ Surface  number 

A1  ■ Radius  of  curvature,  for  R > 6000,  R ■ • 

A2,  A3,  A4  ■ X,  Y,  Z coordinates  of  surface  vertex 

A5  “ Surface  tilt  of  the  vertex  with  respect  to  X-Y  plane 
(in  degrees) 

A6  - Surface  tilt  of  the  vertex  with  respect  to  X-Z  plane 
(in  degrees) 
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b.  N2  < 0 parabolic  or  aconlc  surface 

0 Card  1 

N2  Surface  number 

A1  ■ Radius  of  curvature  (meaningless) 

« A 2,  A3,  A4  ■ X,  Y,  Z coordinates  of  surface  vertex 

A5  ■ Surface  tilt  of  the  vertex  with  respect  in  X-Y  plane 
(in  degrees) 

A6  ■ Surface  tilt  of  the  vertex  with  respect  to  X-Z  plane 
(in  degrees) 

Card  2 

-N2  * Surface  number 

Al,  A2,  A3.  A4,  A3,  and  A6  *■  Aspheric  coefficients 
x - A1  * y2  + A2  ♦ z2  + A3  • r1*  + A4  • r6  + A5  * r8  + A6  • r10 
r2  - y2  + z2 

c.  N2  > 100  ■+•  elliptic  or-  hyperbolic  surface 
N2  - 100  ■ Surface  number 

Al  - RA 

A2 , A3,  A4  ■ X,  Y,  Z coordinates  of  surface  vertex 
A5  - R| 

»•  A6  = Rq 

ix  - Al)j  + jy  - A 2) 2 (2  - A3)2 

• ka  S rg 

R|  and/or  Rq  are  positive  for  an  ellipse  and  negative  for 
a hyperbola . 

• N1  ■ 2 "*■  Read  refractive  Indexes  at  3 wavelengths  and  aperture 
Information  for  each  surface 

M2  ■ Surface  number 

Al  - Refractive  index  for  material  following  surface  N2  at  \1 

A2  ■ Refractive  index  for  material  following  surface  N2  at  X.2 

A3  - Refractive  index  for  material  following  surface  N2  at  A 3 

NOTE:  If  Al,  A2,  A3  are  all  negative,  surface  is  a mirror.  If  N2  ” 

1,  Al  >0,  A2  • 0,  and  A3  - 0,  only  one-color  rays  will  be  traced.  If 
N2  ■ 1,  Al  > 0,  A2  > 0,  and  A3  ■ 0,  only  first  2 colors  will  be  traced. 

A4  • Maximus  aperture  radius  at  this  surface 
A5  - Minimum  aperture  radius  at  this  surface 
A6  ■ Meaningless. 
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Ray  Data  Cards 


Refer  to  Example  3 on  FORTRAN  coding  sheet  (Figure  2) . 

N1  - 3 Read  ray  data  for  individual  rays 

N2  - Ray  number 

A1  - X coordinate  of  ray  N2 

A2  ■ Y coordinate  of  ray  N2 

A3  » Z coordinate  of  ray  N2 

AA  - Angle  ray  makes  with  X axis  in  tangential  plane  in  degrees 
A5  ■ Angle  ray  makes  with  X axis  in  sagittal  plane  in  degrees 

COMPUTATION-TYPE  CARDS 

Refer  to  Example  3 on  FORTRAN  coding  sheet  (Figure  2)  for  format. 

Object  and  Aperture  Information  for  Ray-Trace  Mode 

N1  ■ A Read  object  position  and  generate  ray  data  for  spot  diagram 
. If  0 < N2  < 20  N2  ■ number  of  rings  of  rays  for  snot  diagram 

- tf 2 « -<  N2  > .20  N2  ■ number  of  rays  in  square  lattice 

If  N2  < 0 N2  - number  of  rays  generated  in  one  half  of  the 
i fan  of  rays 

' Al  ■ Range  of  object 

A2  -•  Elevation  In  degrees 
A3  - Azimuth  In  degrees 

1 

AA  * Maximum  aperture  (at  entrance  pupil) 

A5  ■ Minimum  aperture  (at  entrance  pupil) 

A6  “ Stop  position 

Object  and  Aperture  Information  for  Design  Mode 

N1  ■ 5 Read  aperture  'and  , field  angle  and  generate  aberration  rays 
N2  » Nunber  of  rays  • 17  ’ *•  r*  -:r:.  r-. 

Al  - A6  - Identical  to  values  for  N1  ■ A Al  - A6 
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Ray  Trace  Command  Card 

N1  - 6 Compute  ray  trace  from  information  contained  on  N1  - 1,  2, 
and  either  3 or  4 but  not  both. 

N2  - Meaningless 

A1  ■ Increment  of  radius  about  focus.  The  number  of  rays  within 
this  circle  is  contained  in  the  output.  Then  the  radius  is  incremented 
and  rays  are  counted  again,  etc. 

Design  Command  Card 

N1  • 7 Compute  aberrations — automatic  design  mode 
N2  ■ Number  of  variables 

A1  -*■  Defines  number  of  iterations  ■ 1000/A1 
A2  -*■  Defines  Increment  limit  ■■  A2  * DVS 

The  following  two  descriptions  refer  to  special  cards  for  the  J 

design  section. 

Variable  Specification  Card 

One  card  for  each  variable.  Refer  to  Example  4 on  FORTRAN  coding 
sheet  (Figure  2)  for  format. 

DVS  -*■  Used  with  A2  on  N1  - 7 card  to  define  Increment  limit  for 
each  variable  j 

WV  •►1.0  for  all  cases 

KS  - Surface  number 

LS  » Type  of  variable  (what  about  the  surface  the  computer  can 
change) 

LS  r 1 4 radius  of  curvature  ACKg  (C  - 1/r) 

■ 2 ■*  X position  AXjcs 

■ 3,  4,  5 aspheric  constants  A,  B,  and  C,  respectively 

■ 6,  7 aspheric  constants  E and  F,  respectively  (A-F 
eorfcipond  to  2nd,  4th. .. .through  10th  order  terms) 

■ ® allows  a change  In  radius  at  the  same  rate  as  the 
vertex  changes  ► thickness  changes,  power  remains 
the  sane;  Argg  - -AX^ 
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LS  ■ 9 no  thickness  change  distance  between  present  and 
following  surface  is  unchanged;  AXj^g  ■ AX^g+j 

■ 10  changes  shape  of  lens  but  retains  power,  curvatures 
change  at  same  rate;  ACgS  ■ 

■ 11  surfaces  1 and  3 of  a second  surfaced  (Mangln) 
mirror  must  remain  equal;  ACgS  “ ACj^g+2 

MS  - 1 ■+  Meaningless 

m 2 enter  ray  trace — continue  reading 

■ 3 -*■  end  of  design  Input  data,  start  optimizing 


Aberration  Optimization  Selection  Card 


Refer  to  Example  5 on  FORTRAN  coding  sheet  (Figure  2)  for  format. 


FZP  “ Desired  effective  focal  length  of  the  system.  The  following 
variables  correspond  to  factors  that  can  be  optimized.  One  card  specifies 
the  action  taken  on  all  factors.  A factor  set  equal  to  1.0  implies  no 
correction  desired,  higher  values  (i.e.,  10.0)  require  the  program  to 
optimize. 


Factor 

Corrects  for 

WBF 

Back  focal  length 

WFZ 

Effective  focal  length 

WSA 

Spherical  aberration 

WSZ 

Zonal  spherical  aberration 

WPC 

Primary  field  curvature 

WCM 

Coma 

WAT 

Astigmatism 

WST 

Oblique  spherical  (tangential) 

wss 

Oblique  spherical  (sagittal) 

WCL 

Longitudinal  color  (transverse  axial) 

WCT 

Lateral  color 

WZL 

Lower  zonal  spherical 

WZU 

Upper  zonal  spherical 

NWC  TM  2626 


Reticle  Modulation  Calculation  Command  Card 


Refer  to  Example  3 on  FORTRAN  coding  sheet  (Figure  2) . 

N1  ■ 8 Compute  reticle  modulation — static-gain  curves 

N2  ■ 0 -*■  No  data  read  lr- 

N2  - 1 -*■  Read  in  frequency  response  curve  (data  should  be  in 
format  illustrated  in  coding  sheet  Example  7,  Figure  2) 

N2  » 2 -►  Read  in  reticle  pattern  (transmission  in  each  zone, 
(TAU(K))  and  response  curve.  Data  should  be  in  formats 
illustrated  in  coding  sheet  Examples  6 and  7,  respectively. 
Figure  2.  That  is,  the  first  N2  = 2 card  has  the  format  of 
Example  6.  The  next  N2  • 2 card  has  the  format  of  Example  7 
in  Figure  2. 

A1  » Radius  of  reticle 

A2  " length  of  reticle  sector 

A3  ■ Angular  width  of  reticle  sector  (degrees) 

A4  - Incremental  angle  of  rotation  (degrees) 

* 

A5  ■ Rotational  frequency  (cps) 

A6  " Scan  radius 


FIGURE  2.  Sample  of  Fflled-In  FORTRAN  Coding  Sheet 
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DEFINITIONS  OF  ABERRATION  TERMS  USED  IN  THIS  PROGRAM 


Since  different  mathematical  treatments  of  optical  ray  tracing 
give  rise  to  different  meanings  for  the  same  or  similar  aberration  names , 
b brief  definition  for  the  terms  used  in  this  program  follows. 

These  terms  are  explained  using  diagrams  for  on-axis  rays  (Figure  3) 
and  off-axis  rays  (Figure  4),  and  the  definitions  are  grouped  accordingly. 


A 


SYMBOL 


DEFINITION 


"ABERRATION"  NAME 

EFFECTIVE  FOCAL  LENGTH 
PRIMARY  CURVATURE 
COMA 

SECONDARY  CURVATURE 
ASTIGMATISM 

OBLIQUE  TANGENTIAL  SPHERICAL 
OBLIQUE  SAGITTAL  SPHERICAL 
LATERAL  COLOR 

* C AND  F SUBSCRIPTS  REFER  TO 


FIGURE  4. 


EFL 

Ygg/SIN  0 

PC 

*S8-BF1 

CM 

V47  “ VS6 

SC 

xss-,Fi 

AT 

PC  - SC 

OST 

X«7  “ X66 

OSS 

XSS  “ *SB 

TCH 

V8BC  “ VM 

AND  F LIGHT  (DIFFERENT  COLORS!. 


OfT-Axi*  Aberration. 
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METHODS  OF  CALCULATING  ABERRATIONS 


Table  1 describes  the  17  rays  that  are  used  for  calculating  • 

aberrations  In  the  design  mode.  Table  2 lists  the  mathematical 
expressions  used  by  the  program  to  attach  numerical  values  to  the 
system  aberrations.  Numerical  subscripts  in  Table  2 refer  to  ray 
numbers  from  Table  1. 


TABLE  1.  Aberration  Rays. 


Ray  1 

Outer  axial  ray,  \^a 

Ray  10 

502  Zonal  axial  ray,  \^b 

Ray  2 

502  Zonal  axial  ray, 

Ray  11 

Upper  Inner  tangential 
ray,  Xp 

Ray  3 

Inner  axial  ray,  XD 

Ray  12 

Lower  Inner  tangential 
ray,  Xp 

Ray  4 

Upper  outer  tangential 
ray,  XD 

Ray  13 

502  Zonal  axial  ray,  X^c 

Ray  5 

Upper  Inner  tangential 

wy,  \ 

Ray  14 

Upper  inner  tangential 
ray.  Xc 

Ray  6 

Lower  Inner  tangential 

«y,  Ad 

Ray  15 

Lower  Inner  tangential 
ray.  X£ 

Ray  7 

Lower  outer  tangential 
ray,  X,, 

Ray  16 

252  Zonal  axial  ray,  XQ 

Ray  8 

Outer  sagittal  ray,  Xp 

Ray  17 

752  Zonal  axial  ray,  XQ 

Ray  9 

Inner  sagittal  ray,  XQ 

* Xp  ■ Middle  wavelength. 
b Ap  " Short  wavelength. 


Xc  ■ Long  wavelength. 
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TABLE  2.  Mathematical  Expressions  for  Aberrations  (Figure  S). 


Back  focus 


Focal  length 


Spherical  aberration 


Zonal  spherical 
aberration 


Primary  curvature 


Secondary  curvature 


[V 

cos  a~ 
m v ■ 

3 cos 

F 

- Yl/sln  9 

/ cos  a. 

cos  a,  ^ 

l / / cos  a. 

cos  a 

T 

- Iy  4 T, 

o 1 

/ f ^ i 

X 

y 5 cos  5 

6 cos  0^ 

}/  ycos  65 

cos  0( 

Astigmatism 


Oblique  spherical 

aberration 

(tangential) 


Oblique  spherical 

aberration 

(sagittal) 


■> 
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SAMPLE  PROBLEM 


In  this  section  a sample  problem  is  stated  and  the  method  of  solving 
it  is  given.  Then  a sample  output  and  the  interpretation  of  the  output  are 
presented. 


PROBLEM 

Given  a lens  with  the  following  parameters,  find  the  optimum  lens 
shape  to  reduce  coma  while  driving  the  effective  focal  length  to  3.5  inches. 
The  initial  parameters  are: 


I 

R(I) 

XV(I) 

nd(i) 

nf(i) 

NC(I) 

1 

7000 

0.5 

1.0 

1.0 

.1.0 

2 

-1.5 

3.0 

1.51124 

1.5169 

1.50883 

3 

-1.0 

3.35 

1.0 

1.0 

1.0 

A 

7000 

6.801 

1.0 

1.0 

1.0 

METHOD  OF  SOLUTION 

Using  the  design  section  of  342*WCRF .OPTIC,  allow  radii  of  curva- 
ture for  surfaces  2 and  3 to  be  variables  while  choosing  to  correct  for 
^ EFL  and  coma. 

The  data  cards  required  to  perform  the  above  calculations  are 
illustrated  in' Figure  6. 

INTERPRETATION  OF  OUTPUT 

A sample  output  data  printout  consisting  of  four  sheets  is  shown 
in  Figure  7.  Sheet  1 displays  the  input  data  which  describes  the  initial 
system  plus  the  output  for  each  iteration  of  the  design  mode  section.  The 
tint  group  of  columns  represents 

I “ Radius  of  curvature 
XV  m Position  of  vertex  along  x (optical)  axis 
TV  - Position  of  vertex  along  y axis  (decentering  direction) 

ZV  - Position  of  vertex  along  s axis  (decentering  z direction) 

Gamma  T - Surface  tilt  in  tangential  (x-y)  plane 

Gamma  S - Surface  tilt  in  sagittal  (x-z)  plane 
ND  - Refractive  index  in  D light 
( MF  ■ Refractive  index  in  F light 

NC  - Refractive  index  in  C light 
Surface  number 


23 


FIGURE  6.  Sample  Input  Data  Coding  Sheet 
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The  second  grouping  of  columns  on  sheet  1 refers  to 

AY,  AZ  ■*  Second-order  aspheric  terms  (x-y  and  x-z  planes, 
respectively) 

B,  C,  D,  E ■ 4th-,  6th-,  8th-,  and  lOth-order  aspheric  terms, 
respectively 

H Max.  “ Maximum  aperture  at  each  surface  (measured  from 
the  vertex  to  outer  edge  of  the  element) 

H Min.  ■ Minimum  aperture  at  each  surface  (measured  from 
the  vertex  to  inner  edge -of  the  element) 

The  following  information  appearing  on  6heet  1 of  Figure  7 is 
output  from  the  design  section.  Refer  to  Table  2 for  definitions  of 
these  terms. 

BF  ■ Indicates  correction  for  back  focus 
EFL  ■ Effective  focal  length 
SA  - Spherical  aberration 
ZSA  - Zonal  spherical  aberration 
PC  ■ Primary  curvature 
CM  ■ Coma 
AT  ■ Astigmatism 

OST  ■ Oblique  spherical  (tangential) 

OSS  - Oblique  spherical  (sagittal) 

LCH  “ Longitudinal  (transverse  axial  color) 

TCH  “ Lateral  color 

ZSL  ■ Lower  zonal  spherical 

ZSU  ■ Upper  zonal  spherical 


t 
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The  numbers  directly  under  these  columns  are  the  values  of  these 
terms  at  each  iteration.  The  numbers  under  each  of  these  rows  and 
Indented  somewhat  are  the  values  the  program  generated  for  these 
variables,  the  .user  specified  as  changeable. 

In  general.  Iterations  will  continue  until  the  user  specified 
limit  la  reached  unless  a solution  for  the  problem  Is  reached  beforehand. 
Occasionally  though,  the  terms  chosen  to  be  minimized  can  work  In 
opposition — l.e.,  as  one  Increases,  the  other  reduces.  This  situation 
can  cause  an  oscillation  in  the  values  of  the  aberration  terms.  When 
this  happens  it  Is  quite  possible  that  the  best  solution  Is  not  the 
l*st  Iteration  results,  Instead  an  earlier  Iteration  may  display  the 
optlsuan  configuration. 
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In  Che  last  three  rows  of  the  design  output  one  finds  the  final 
Iteration  term  values  in  the  first  row  and  numbers  which  are  related  to 
the  rates  of  change  in  these  term  values  in  row  two.  Zn  row  three  are 
numbers  that  refer  to  the  relative  weighting  of  aberrations  to  be 
corrected. 

Sheet  2 of  the  sample  output  (Figure  7)  and  subsequent  sheets  are 
the  output  of  the  ray-trace  mode.  The  first  two  groupings  of  columns  on 
sheet  2 contain  the  data  that  was  input  to  the  ray-trace  section.  In 
this  example,  the  design  section,  which  always  automatically  calls  for  a 
ray  trace,  uses  data  calculated  in  the  last  design  iteration  for  input 
ray-trace  data.  (Note  the  change  in  R2  and  made  by  the  design  section.) 

The  next  grouping  of  columns  contain  data  for  each  ray  at  each 
surface.  A diagram  illustrating  these  terms  and  definitions  of  the 
terms  is  presented  in  Figure  8. 

The  unlabeled  three-digit  column  of  integers  defines  the  ray  number 
and  surface  number  for  each  row.  For  example,  101  means  the  row  contains 
data  for  the  first  ray  and  the  first  surface,  904  means  the  row  contains 
data  for  the  ninth  ray  and  the  fourth  surface.  The  two  remaining  unlabeled 
columns  represent  back  focus  and  effective  focal  length. 

When  the  design  section  is  used,  the  output  ray  data  is  printed  for 
nine  rays  at  each  surface  of  the  system.  After  the  ray  data  is  printed, 
three  rows  of  data  with  eight  columns  in  each  row  are  printed.  Row  one 
consists  of: 

1.  Range  of  object 

2.  Elevation  of  object  in  degrees  below  the  horizontal  (x-y  plane) 

3.  Aximuth  angle  in  degrees  from  the  horizontal  (x-z  plane) 

4.  5.  and  6.  The  x,  y,  and  z positions  of  the  energy  centroid 

at  the  image  plane  (meaningful  only  for  individual  ray  traces 
not  for  design  called  ray  traces) 

7.  Nimber  of  rays  reaching  last  surface  (check  for  vignetting 
and  internal  reflections) 

. Diameter  of  the  root  mean  square  spot  size 
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Row  two  and  row  three  are  related.  Each  column  in  row  two  repre- 

* sents  a radius  about  the  energy  centroid;  in  row  three,  below  each  radius 
is  a number  which  describes  how  many  rays  are  contained  within  the  annulus 
formed  by  the  previous  column's  radius  and  the  present  column's  radius. 

* In  the  first  column  case  the  number  of  rays  refer  to  the  number  within 
the  circle  defined  by  the  first  column's  radius  about  the  centroid.  Rows 
two  and  three  are  not  meaningful  when  the  design  section  initiates  a rav 
trace. 

Below  the  three  rows  of  data  Just  explained  is  a second  grouping  of 
ray-trace  data.  This  output  data  contains  the  results  of  the  ray-trace 
mode  when  the  second  refractive  index  specified  on  the  NS  - 2 card  is  used. 
In  other  words,  the  data  is  information  about  rays  traced  through  the 
system  for  a different  color. 

The  three  rows  following  this  group  of  ray-trace  data  are  the  same 
as  the  first  group  of  three  rows,  with  the  exception  that  the  numbers  of 
rays  contained  within  the  annuluses  are  cumulative  for  the  two  colors 
I The  same  is  true  for  the  three  rows  following  the  ray- trace  data  for  the 

third-color  rays. 
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FIGURE  7.  Simple  Output  Date  (Sheet  3 of  4). 
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H - DISTANCE  OF  RAY'S  INTERSECTION  WITH  SURFACE  FROM  VERTEX. 

X • POSITION  ALONG  X-AXIS  OF  A PERPENDICULAR  DROPPED  FROM  THE  INTERSECTION 

OF  THE  SURFACE  AND  THE  RAY.  : 

V * RAY  HEIGHT  MEASURED  ALONG  THE  ABOVE-OESCRIBED  PERPENDICULAR  IN  Y-PLANE. 

Z * RAY  HEIGHT  MEASURED  ALONG  THE  ABOVE-OESCRIBED  PERPENDICULAR  IN  Z-PLANE. 

THETA  T - RAY  ANGLE  AFTER  INTERSECTION  WITH  THE  SURFACE.  MEASURED  IN  DEGREES  OFF 
THE  HORIZONTAL  IN  TANGENTIAL  PLANE. 

THETA  S - SAME  AS  THETA  T.  BUT  IN  SAGITTAL  PLANE.  ' J 

i ? 
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FIGURE  8.  Illustration  of  Ray-Data  Terms. 
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TROUBLESHOOTING 


Occasionally  attempts  to  solve  new  problems  will  result  in  error 

* messages  whose  meanings  are  not  immediately  obvious.  The  following  list 
of  errors  and  corrections  is  presented  to  give  a few  starting  points  in 
problem  debugging.  At  the  end  of  the  list  is  Included  one  problem  that 
does  not  produce  an  error  message. 

Error:  "No  plotter  initialization  routine  has  been  called .. ." 

Check:  N1  - 6 or  N1  ■ 7 cards.  Hake  sure  either  one  or  the  other 
la  present  and  that  N1  is  contained  in  column  4. 

Error:  "Ray  height  exceeds  maximum  aperture  at  interface  No.  ... 

Check:  Y values  versus  H values  at  each  X location  for  the  ray 
in  question  as  follows: 

1.  If  Y and  H values  are  nearly  equal  and  the  X value 
la  near  the  specified  surface  location,  a simple 
re-evaluation  of  the  surface  aperture  size  (N1  *■  2 
card),  field  angle,  entrance  pupil  diameter,  or 

* limiting  aperture  location  (all  specified  on  either 

the  HI  - 4 or  5 card)  may  be  required. 

* 2.  If  Y and  H values  differ  significantly  and  if  the  X 

value  deviates  from  the  specified  Interface  location, 
check  the  following: 

a.  The  surface  cards  and  numbering  to  see  that  they 
are  ordered  in  the  sequence  that  the  rays  are 
expected  to  reach  the  surfaces.  Watch  out  for 
negative  lens  thicknesses  as  a result  of  design 
mode  changes. 

b.  The  possibility  that  a ray  may  Intersect  the 
surface  on  the  wrong  side  of  the  curve  as  illus- 
trated in  Figure  9(a)  and  (b) . To  prevent  this. 
Insert  a dummy  plane  surface  immediately  before 
the  surface  exhibiting  this  difficulty. 

Error:  "**  Points  out  of  Range  **" 

Check:  Craphlcs  scale  limits  specified  on  the  card  preceding 
HI  ■ 1 cards.  Points  can  be  out  of  range  in  either  the 
X or  Y direction.  Also  check  which  surface  was  specified 
for  the  starting  point  of  the  graphics  output. 


MM 
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<•>  ERRONEOUSLY  REFRACTED  RAY 


M SITUATION  CORRECTED  BY  INTRODUCING  DUMMY  SURFACE  NO.  7 


FIGURE  9.  Illustrations  of  Ray  Intsmctioa. 
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Error:  "Ray  internally  reflected  at  Interface  ..." 

Check:  For  tight  curvatures  and  high  Indexes  of  refraction. 

Adjust  shapes  or  focal  lengths  to  accommodate  rays . 

Error:  "The  interpretation  of  meaningless  data  was  attempted." 

Check:  Card  formats,  card  sequencing,,  characters  punched  on  cards. 

Error:  "Ray  does  not  intersect  interface  number  ..." 

Check:  Position  of  sequential  elements  along  optic  axis.  Rays 

will  not  reverse  directions  unless  a mirror  reflects  them. 
Curvature  too  tight. 

Problem:  Terms  to  be  minimized  are  not  improving  or  are  not  improving 
fast  enough. 

Check:  1.  Formatting  on  aberration  minimization  specification  card. 

2.  Adequate  number  of  variables  (one  for  each  aberration  to 
be  minimized  and  always  2 2) . 


3.  Large  enough  Increment  limits  for  each  variable. 

4.  Nunber  of  possible  Iterations. 

5.  Variables  selected  are  consistent  with  aberrations  to 
be  minimized.  For  instance  to  correct  effective  focal 
length  one  would  not  choose  to  change  lens  shape  while 
retaining  the  same  power. 
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